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In the absence of a confining potential, the boson Hubbard model in its ground state is known 
to exhibit a superfluid to Mott insulator quantum phase transition at commensurate fillings and 
strong on-site repulsion. In this paper, we use quantum Monte Carlo simulations to study the 
ground state of the one dimensional bosonic Hubbard model in a trap. We show that some, but not 
all, aspects of the Mott insulating phase persist when a confining potential is present. The Mott 
behavior is present for a continuous range of incommensurate fillings, a very different situation 
from the unconfined case. Furthermore the establishment of the Mott phase does not proceed via a 
quantum phase transition in the traditional sense. These observations have important implications 
for the interpretation of experimental results for atoms trapped on optical lattices. Initial results 
show that, qualitatively, the same results persist in higher dimensions. 
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A considerable amount of work has been done in the 
last decade to determine the ground state phase diagram 
of correlated bosons on a lattice described by the "boson- 
Hubbard" Hamiltonian On-site repulsion can pro- 

duce a Mott insulating phase at commensurate fillings, 
with a quantum phase transition to a superfluid as the 
density is shifted or the interaction strength weakened. 
Longer range interactions can cause charge density wave, 
stripe, or even supersolid order Extensions to dis- 

ordered systems have allowed the detailed study of the 
interplay of randomness and interactions in quantum sys- 
tems 1). 

Recently, the trapping of atoms on optical lattices has 
given another experimental reahzation of these bosonic 
phases. However, the quadratic confining potential 
present in addition to the regular "lattice" potential, 
leads to a number of fundamentally new, and open, is- 
sues: (i) Does the confining potential preclude the for- 
mation of Mott regions by providing a continuous, un- 
bounded, distribution of local site energies? (ii) If an 
insulating phase still exists, how is it characterized? 
(iii) What are the quantitative values of the trap curva- 
ture and interaction strength that support Mott phases? 
These questions are largely unaddressed in the literature. 

In this paper we report the first quantum Monte Carlo 
simulation of the one-dimensional boson-Hubbard model 
in a confining quadratic potential and provide a quan- 



titative map of the state diagram. To our knowledge, 
the only non-mean-field work on this problem [|ll[ is on 
very small systems (five particles) . We find that the trap 
changes the physics fundamentally from that found in 
earlier simulations and subsequent analytic |^ studies. 
For example, the vanishing of the global compressibility, 
discussed at length in and reported in recent mean- 
field studies of the boson-Hubbard model in the context 
of optical lattices , but which ignores the confining po- 
tential, is absent. Other recent papers ^-[1^] on bosons in 
optical traps have likewise emphasized similarities to the 
physics in the absence of a confining potential, and have 
used values of the unconfined lattice critical coupling to 
compare with experimental data. 

We find the inhomogeneous potential resulting from 
the confining trap is a crucial feature in discussing Mott 
regions at incommensurate fillings j|] , and the physics of 
bosons on optical lattices generally. The Mott insulat- 
ing regions exist above a threshold interaction strength, 
even without the commensurate filling required in the 
non-confined case. This is a consequence of the inho- 
mogeneous distribution of boson densities which allows 
extended Mott domains with, for example, one parti- 
cle per site to coexist with regions of other density and 
is a unique feature which distinguishes the behavior of 
the confined model. Similarly, the global compressibil- 
ity is non-zero for all densities, including commensurate 
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ones. The different regimes of the state diagram can be 
characterized by the topology of the density and local 
compressibility profiles. Therefore, in this context, the 
(global) density, p, defined as the total number of parti- 
cles divided by the system size, loses its meaning since 
the particles are not uniformly distributed and the "sys- 
tem size" itself is ill-defined. Adding particles can push 
bosons deeper into the confining potential at the edges 
of the system thus changing the "size" . 

We review briefly the properties of the ground state of 
the d-dimensional non-confined boson Hubbard model. 
This model has two phases: a Mott insulator at commen- 
surate fillings and sufficiently strong interactions, and a 
superfluid elsewhere . The critical behavior is of two 
types: mean field for transitions induced by tuning the 
density, and of the (d -I- 1) dimensional XY universal- 
ity class when the interaction strength is swept at fixed 
commensurate filling. One of the key new results of this 
paper will be that this special status of commensurate 
filling is lost in the case of a confining potential since 
commensuration itself well defined only locally.. 

We study the Hamiltonian 

H = -tJ2iaUi+i + at+i^t) (1) 

i 

+ VoY^ n,(n, - f ) + 14 Y^(i - Ljlf n„ 

i i 

at zero temperature. Here t measures the boson kinetic 
energy, Vb the on-site repulsion, Vc the curvature of the 
quadratic confining potential and L the number of sites. 
In the presence of the trap, the value of L should be cho- 
sen such that for the given trap curvature, the bosons do 
not see the edge of the system and therefore do not leak 
out. Our simulations were done with the world-line quan- 
tum Monte Carlo algorithm in the canonical ensemble 
|^,|^. The chemical potential /i — dE/dN is obtained 
by differentiating numerically the energy with respect to 
the particle number [T^ ]. In the presence of a confining 
potential, it is important to measure the local density of 
bosons, rii — {a\ai) as well as the local compressibility. 

Figure |l| shows the evolution of the local boson den- 
sity with increasing total occupancy of the lattice. At 
low fillings the density profile is smooth, with an inverted 
parabolic shape reflecting the confining potential. Above 
a critical filling of about 30 bosons (for this choice of 
Vb and Vc) a plateau with a local filling of one boson 
per site develops in the density profile which is analo- 
gous to the Mott structure of N vs fj. in the unconfined 
model. This plateau indicates the presence of an incom- 
pressible, insulating region where defined above, 
drops to a small but finite value (Fig. ||) which vanishes 
for Vb ^ oo. Here, this behavior of Ki will be taken as the 
signal for a Mott region. As the density is increased fur- 
ther, the plateau widens spatially. But when the energy 




FIG. 1. The evolution of the local density m as a function 
of increasing total number of bosons. The trap curvature is 
Vc = 0.008, L = 100 and the on-site repulsion is Vb = 4. At 
low fillings the system is in a superfluid phase. Mott insulat- 
ing behavior appears as the density is increased, but then at 
yet larger fillings a superfiuid begins to form at the center of 
the insulating region. 

cost of extending the plateau to increasingly large values 
of the confining potential becomes prohibitive, the occu- 
pancy begins to exceed one at the center of the lattice, 
indicating a breakdown of Mott behavior there, but not 
everywhere. Increasing the filling further, for example 
Ni, = 116, eventually produces a second Mott region in 
the center of the system with a local filling of two bosons 
per site without destroying totally the first Mott region. 
Four slices from Fig. |l] are shown in Fig. |^ along with the 
local compressibilities. It is clear that at higher boson 
numbers, richer structures where the local compressibil- 
ity vanishes at several locally commensurate densities can 
occur. Mean-field work in two dimensions ||ll| shows a 
similar coexistence of Mott and superfluid regions. 

A central feature of the Mott phase transition of the 
unconfined boson-Hubbard model is global incompress- 
ibility: A charge gap opens up, i.e. the density gets 
"stuck" at p = integer for a range of chemical potentials 
fi. One might, then, crudely interpret the spatial depen- 
dence of the local density in the confined case. Fig. H, as 
rather analogous to the chemical potential dependence 
in the unconfined case |pd| . This assumption, and the p 
versus fi curve in the unconfined case, allow us to calcu- 
late the site at which a Mott domain is entered or exited. 
These are shown as arrows in Fig. ^. However it is vital to 
emphasize that while the confined system has locally in- 
compressible regions, the global compressibility is never 
zero which is seen clearly in Fig. ^ The main figure 
should be contrasted with the nonconfined case (inset). 

An important difference in the behavior of the local Ki 
is especially evident in one dimension where, in the un- 
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FIG. 2. Cuts across Fig. 1 show the compressibihty profile 
Hi (sohd line) associated with the local density rii (circles). 
The fillings are Nt — 25(a), 33(a), 50(a) and 60(a). Ki is very 
small when rii — 1. For the arrows see text. 
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FIG. 3. Nb (number of bosons) as a function of chemical 
potential, /i for Vo = 4.5. No globally incompressible Mott 
plateau is observed. Inset shows the unconfined case. 

confined case, the global compressibility diverges as 
the Mott lobe is approached, k oc \p — for the first 
lobe. Here, instead, we find Hi oc {rii ~ 1), as shown in 
Fig. 1^. The origin of these differences is of course that the 
global compressibility, k = l3{{J2ij niUj) — (n)'^), probes 
density correlations at all length scales. In the unconfined 
case, contrary to the confined system, the establishcment 
of the Mott phase is a true quantum phase transition: It 
happens collectively throughout the system and the cor- 
relation length diverges. There is, however, an interesting 
"universality" in the trapped system as a Mott region is 
approached. That is, the values of the same even 

as the total filling and the on-site repulsion are varied 
(see Fig. ^). The same behavior is observed for the n — 2 
locally incompressible phase. 

Sets of runs such as those shown in Fig. 1 allow us to 
determine the state diagram as a function of boson filling 
and interaction strength for a given trap curvature. This 
is shown in Fig. ||. Because of the absence of true phase 
transitions, we have referred to Fig. ^ as a state diagram 
rather than phase diagram. As the filling is increased 
at fixed interaction strength, one crosses from a smooth 
density profile to one which has locally incompressible 
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FIG. 4. The local compressibility as a function of local 
density for Vo = 4 and Nh = 35 (o) and A^i, = 80 (o) and 
for Vo = 4.5, Nb = 90 (A) and Nb = 141 (y). « decreases 
linearly with n as the Mott region, n = 1, is approached from 
below or above. This behavior is dramatically different from 
the unconfined system where k diverges. 
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FIG. 5. The state diagram of correlated bosons in a 
quadratic confining potential. The solid lines are to guide the 
eye, the dashed lines are extrapolations. See text for details. 



"Mott" domains, if Vq is large enough. Further increase 
in the filling ultimately leads to the formation of regions 
at the well center where > 1. In Fig. ||, region A ad- 
mits only n & locally incompressible regions with rii = 1, 
as in Fig. g(b) and (c). Region B has > 1 surrounded 
by incompressible regions. Fig. ||(d). Region C is where 
the central part of the system has an = 2 incompress- 
ible region which, when the boundary of the system is 
approached, falls off to a shoulder of = 1 "Mott" re- 
gion before reaching zero density. The Ui = I and Ui = 2 
incompressible regions are separated by compressible re- 
gions. Region D is where the center of the system is 
compressible > 2, bounded by = 2 which in turn is 
bounded by = 1 incompressible regions. Region E has 
no incompressible regions. 

Note that the values of Vq at which the A and C re- 
gions in Fig. ^ are entered are of the same order as those 
of the first two true Mott lobes in the nonconfined case 
1^. This is consistent with the experimental results on 
three-dimensional optical lattices |^ which appear to be 
in agreement with the expected value in the nonconfined 



3 



case. Furthermore, the narrowness of region C could help 
understand why the experiments have not shown signs 
of the n = 2 Mott region, even though rii ~ 2.5 in the 
core of the system. 

One of the interesting experimental results is how 
rapidly coherence is re-established when Vb is suddenly 
reduced from a value large enough to have produced large 
incompressible regions. It was argued that the character- 
istic time is of the order of the tunneling (hopping) time 
between sites |^] . This is entirely consistent with the pic- 
ture we present here. If, for example, the density profile 
is as in Fig. ||(c) when Vq is suddenly reduced to a very 
small value, the system will evolve to a profile like in||(a) 
albeit with higher local density in the center. This is ac- 
complished by particles near the edges hopping towards 
the center, and is greatly accelerated by the fact that 
the trap is much lower near the center than near edges. 
In addition, since there is only one center of nuclcation 
(the geometric center of the system) there is no slowing 
down due to competition at domain walls where different 
nucleation zones meet. 

In summary, we have discussed the nature of lo- 
cally incompressible "Mott" insulating behavior in a one- 
dimensional system of interacting bosons in a confining 
potential at T = 0. We conclude that because of the 
destruction of translation invariance, great care should 
be taken in drawing on the analogy with the unconfined 
case at a fundamental level. For one thing, it can sup- 
port "Mott" behavior off commensurate fillings. While 
"Mott" regions still exist, the critical properties are com- 
pletely altered: Incompressible regions are established in 
a very localized way and not at all critically in the usual 
sense. These localized regions grow or shrink with Vq, 
but are always in co-existence with other regions, some 
compressible and some which might be incompressible 
but at higher local integer filling. In that sense, the for- 
mation of these "Mott" regions is not a true quantum 
critical phenomenon as it is in the unconfined case. 

While we have focussed on the new qualitative physics 
which results from the confining potential, it is impor- 
tant to emphasize that experiments on one and two 
dimensional trapped systems are currently underway 
^ , p7| -[l9|. For these, our paper should provide specific 
quantitative predictions for the critical ratios of interac- 
tion strength to kinetic energy and trap curvature, as a 
function of density. We are currently undertaking these 
comparisons. In addition, in the absence of traps, the 
phase diagram is qualitatively the same in one and two 
dimensions . We expect this to be true in the con- 

fined case too. In fact, initial simulation results in two 
dimensions show this to be true. 

Note Added: A QMC study of the three-dimensional 
system appeared recently |2^] in which the authors also 
conclude, like we do, that one cannot characterize glob- 
ally the transitions discussed here. They discuss a signal 
that can be used experimentally to study the transition. 
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